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: ABSTRACT 

a 

SN 2003dh, one of the most luminous supernovae ever recorded, and the one 
^ ■ with the highest measured velocities, accompanied gamma-ray burst 030329. Its 

■ rapid rise to maximum and equally rapid decline pose problems for any spheri- 
cally symmetric model. We model the supernova here as a very energetic, polar 
explosion that left the equatorial portions of the star almost intact. The total 
progenitor mass was much greater than the mass of high-velocity ejecta, and 
the total mass of ^^Ni synthesized was about 0.5 solar masses. Such asymme- 
tries and nickel masses are expected in the collapsar model. A "composite two- 
dimensional" model is calculated that agrees well with the characteristics of the 
observed light curve. The mass of ^^Ni required for this light curve is 0.55 Mq 
and the total explosion energy, 26 x 10^^ erg. 
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1. INTRODUCTION 

On March 29, 2003 one of the brightest gamma-ray bursts (GRBs) in history was dis- 
covered and locahzed by the HETE-2 sateUite. Astronomers worldwide watched and within 
days were rewarded by the discovery of a Type Ic supernova, SN 2003dh, in precisely the 
same location (Stanek et al. 2003; Hjorth et al. 2003). From the light curve and spectrum, a 
temporal coincidence with the burst was also estimated to be < 2 days. It is now generally 
agreed that the supernova and the GRB came from the same explosion. 

The V-band light curve and spectrum of SN 2003dh closely resembled that of another 
famous supernova-GRB pair, SN 1998bw (Galama et al. 1998) and GRB 980425, but with 
several important distinctions: 1) GRB 980425 had an observed equivalent isotropic energy 
in gamma-rays roughly four orders of magnitude less than GRB 030329; 2) SN 2003dh rose 
to maximum in less than 10 days, SN 1998bw took 16 days (Woosley, Eastman, & Schmidt 
1999); 3) SN 2003dh exhibited higher expansion speeds, up to 40, 000 km s~^; and 4) the 
total energy of ejected relativistic matter in SN 1998bw was < 3x10^° with v > 0.5 c (Li & 
Chevalier 1999). This is less even than the energy in gamma- rays from GRB 030329, unless 
the beaming angle is very small. 

It is generally agreed that SN 1998bw was a very asymmetric explosion whose high 
velocities may not have characterized ejecta at all angles (e.g., Mazzali et al. 2001; Maeda 
et al. 2003), and the same seems hkely to be true of SN 2003dh. Despite the fact that its 
radiation is not beamed like a GRB, an asymmetric supernova can have a quite different light 
curve from a spherical one of the same total energy. The higher velocities in one direction 
lead to the earlier escape of radiation and, in the case of a radioactive power source, a rapid 
decline in the deposition efficiency of gamma-rays. If the radioactivity itself is mixed out 
preferentially along one axis, the efficiency of gamma-deposition is also affected, leading to 
a more rapid rise and decline of the luminosity. On the other hand, the less rapidly moving 
matter ejected in other directions, can continue to contribute an extended tail on the light 
curve, powered by the remaining radioactivity. 

In a pioneering study of light curves from asymmetric explosions, Hoflich, Wheeler, 
& Wang (1999) modeled the light curve of SN 1998bw, but attempted to fit it into the 
general family of Type Ib/c supernovae. Their parameters were thus typical of these common 
events: ^^Ni mass (0.07 to 0.2 Mq), ejected mass (2 Mq), and explosion energy (2 xlO^^ erg). 
They, and Woosley, Eastman, & Schmidt (1999), championed the idea that asymmetrically 
exploded supernovae of the same energy may have different light curves. However, in light 
of SN 2003dh, neither went far enough. 

Here we model the supernova as essentially two components - a slowly moving, high-mass 
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equatorial ejection and a lower mass polar ejection (see also Maeda et al. 2002, 2003). The 
juncture between these is almost discontinuous. A smooth transition would give too broad 
a light curve. The total mass of the progenitor, ~ 10 Mq of helium and heavy elements, the 
^^Ni mass synthesized, ~ 0.5 M©, and the total explosion energy, ~ 10^^ erg, are all quite 
atypical for ordinary Type Ib/c supernovae. 

2. EXPLOSION MODELS 

A series of explosions was calculated for two Wolf-Rayet (WR) stars of final mass 8.39 Mq 
and 15 Mq. Both models had an initial mass of 15 Mq and an initial composition appropriate 
to the helium core of a star with 0.1 solar metallicity. Both models were started with a surface 
rotation rate corresponding to 30 % Keplerian at the equator. It was assumed, however, that 
the lower mass star was in a binary system and lost its hydrogen envelope to the companion 
star during the expansion phase after core hydrogen depletion (Case B mass transfer). This 
model continued to lose mass as a WR-star at a rate given by Braun (1997), reduced by 
^/[Fe/FeQ] (i.e., about a factor 3), and another factor 3 to account for clumping (Hamann 
& Koesterke 1998). The higher mass model, while lacking a hydrogen envelope, assumed no 
further mass loss. This was an artificial way to create two rapidly rotating WR-stars with 
a range of masses. For the given starting mass, the one evolved with mass loss is probably 
the more realistic. 

In both calculations the effects of rotationally induced mixing were included (Heger, 
Langer, & Woosley 2000) throughout the evolution and both reached iron core collapse with 
sufficient angular momentum to form a Kerr black hole and an accretion disk. The adopted 
reaction rate for the ^^C(a,7)^^0 reaction was 1.2 times that of Buchmann (1996). Both 
models had very fine surface zoning, down to less than 10^^ g. The radii of the two stars at 
the time their cores collapsed were 8.1x10^° cm (with mass loss) and 8.8x10^° cm (without 
mass loss). Some other properties of the presupernova stars are given in Tables 1 and 2. 

Explosions were simulated by placing a piston at the location of a large entropy jump, 
around S/NAkB = 4, in the presupernova star when its peak infall velocity had reached about 
1000 km s~^. Such a large change in entropy typically corresponds to a sudden decrease in 
density at the base of the oxygen-burning shell where mass bifurcations often develop in 
supernova models. The star outside this piston was first allowed to collapse to 500 km at 
one-fourth the free fall acceleration (inward movement of the piston). For the lower mass 
model, which only formed a small iron core (Table 2), the entropy jump and piston were at 
1.462 Mq. For the higher mass model, a large collapsing low-entropy core formed (Table 2) 
and two piston locations were explored. One was located at 1.93 Mq, at the edge of the iron 
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core, the other at 2.75 M© where the large entropy discontinuity occurred. 

After the minimum radius was reached, the piston was moved outward supersonically, 
decelerating at constant fraction of gravitational acceleration until it coasted to a halt at 
10,000 km. The initial velocity and deceleration of the piston were adjusted to give the 
desired explosion energy. (The right piston acceleration has been determined using a modified 
regular FALSI algorithm.) This explosion energy here is defined as the final kinetic energy 
of the ejecta for an explosion into vacuum. 

A wide range of kinetic energies was explored for these isotropic explosion (Table 3), 
including energies all the way up to 1.6 x 10^'^ erg. This energy, half the binding energy 
of a typical neutron star, is far more than expected from any realistic neutrino-powered 
model, but, as we shall see, might be relevant for a small amount of mass ejected in a very 
asymmetric explosion in which the energy source is not neutrinos, but gravitational energy 
from accretion into a black hole. Nucleosynthesis was calculated as described in Weaver, 
Zimmerman, & Woosley (1978). An interesting result was the observation of a maximum 
mass of ^^Ni, regardless of explosion energy and depending only on the mass of the progenitor 
and depth of the piston. 

A near constant mass, 0.2 ± 0.05 Mq was synthesized in the 8.39 models (Table 3). 
The near constancy of this upper bound is a consequence of the fast expansion and increasing 
entropy in the center of the ejecta that leads to freeze-out of ^He rather then ^^Ni (Fig. 1). 
Eventually, turning up the energy only increases the ejection of a-particles, not of ^^Ni. 
Because the location of the piston cannot be much deeper in the 8.39 Mq model and because 
any other way of exploding the star, say by energy deposition rather than a piston, would 
give more photodisintegration, ~ 0.2 M0 is the maximum ^^Ni that can be synthesized in a 
spherically symmetric explosion of this star. A larger amount can be made in the 15 Mq star 
without mass loss because more mass sits closer to the piston. For a piston mass of 2.75 Mq, 
the limiting mass of ^^Ni was ~ 0.6 Mq. For a piston mass of 1.93 Mq, the edge of the iron 
core, the limiting ^^Ni mass was 1.1 Mq (Table 3). One might get still larger ^^Ni masses by 
going to helium stars above 15 Mq. However, this is the largest helium core one expects for 
stars near solar metallicity and larger mass cores will require even greater energies to expand 
rapidly enough to explain SN 2003dh. 

That these limits are within a factor of two of the mass of ^^Ni inferred (see below) for 
both SN 1998bw and SN 2003dh is interesting, but probably coincidental. The existence of 
an upper bound for shock powered models has interesting implications though (§ 5). 
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3. MODEL LIGHT CURVES 

Light curves were calculated for all the explosion models using the Kepler code as 
described, e.g., in Woosley, Eastman, & Schmidt (1999). Gamma-rays from ^^Ni and ^^Co 
decay were assumed to deposit locally. The gamma-ray opacity was 0.037 cm^g~^. For 
the diffusing radiation, opacity was assumed to be predominantly electron scattering. The 
electron density was calculated by solving the Saha equation for the ionization structure at 
each point. Though calculated assuming a single-temperature, flux-limited diffusion, and 
a simple model for gamma-ray deposition, the curves should be qualitatively correct and 
suffice for present purposes. 

Fig. 2 and the first frame of Fig. 3 show the light curves expected when the supernova 
experiences "moderate" mixing. Mixing was simulated by a running average of the composi- 
tion within a region of mass. Am, set here to 10 % of the final mass of the star. That is, the 
actual composition gradients were smoothed by averaging within a defined band of masses, 
and this band was itself moved out, zone by zone, from the piston to the surface. This pro- 
cedure was repeated four times in each model. For the 8.38 Mq model with 1.25 x 10^^ erg 
explosion energy, for example, this resulted in ^^Ni being mixed far enough out that it had 
50% of its central value (i.e., just above the piston) at about 3M0 and 10% at 5.5 M© 
(as measured from the center of the collapsed remnant). The second panel of Fig. 3 shows 
similar light curves for the 8.35 Mq model when the composition is completely homogenized, 
i.e., made to be the same from center top surface. This might be the case if a vigorous 
asymmetric flow is responsible for exploding the star. 

These are bolometric light curves, which properly should be compared only with lumi- 
nosities integrated across the UVOIR bands. However, based upon our experience with SN 
1998bw, we shall compare them with the observed V-band light curve of SN 2003dh (Hjorth 
et al. 2003). This shows the brightness rising to maximum at 10 - 13 days (rest frame) and 
declining by 1 magnitude when the supernova was about 30 days old. The authors further 
estimate that the supernova was "slightly brighter" than SN 1998bw for which a peak lumi- 
nosity 1.0x10^^3 erg s"^ (Woosley, Eastman, & Schmidt 1999) and ^^Ni mass 0.3 to 0.4 M© 
has been estimated (Nakamura et al. 2001; SoUerman et al. 2002). 

Though the ^^Ni mass ejected varies with kinetic energy and progenitor mass (Table 3), 
the supernova must have made a single value. Knowing the approximate luminosity we 
needed at peak, another series of light curves was calculated for which the yield of ^^Ni 
was normalized to 0.5 M© in all the explosions, even for the low mass progenitor. This is a 
typical value for the very energetic explosions (Table 2), is close to the 0.3 Mq (Sollerman et 
al. 2002) to 0.4 Mq (Nakamura et al. 2001) of ^"^Ni inferred for SN 1998bw, and is consistent 
with the observation that SN 2003dh was a little brighter than SN 1998bw (Hjorth et al. 



- 6 - 



2003). The fact that it is allowed to exceed the hydro dynamical limit derived in the previous 
section is justified because of the alternative way ^^Ni is made in the coUapsar model (§ 4). 

Such energetic models with so much ^^Ni are likely to be thoroughly mixed, more so 
than usual supernovae. Mixing proved necessary to obtain a good fit to the light curve of 
SN 1998bw, which accompanied GRB 980425 (Chugai 2000), and is favored by the rapid rise 
time in SN 2003dh. The third panel in Fig. 3 shows the light curves expected for the 8.39 Mq 
models when the ejecta are thoroughly mixed and forced to contain a fiducial 0.5 Mq of ^^Ni. 
This was accomplished by replacing as much of the material directly above the piston by 
pure ^^Ni as was needed to obtain the desired total ^^Ni mass. The replacement was done 
100 s after the explosion when essentially all thermonuclear reactions have ceased, but before 
any mixing was applied. 

Fig. 4 shows the terminal velocity for both models as a function of mass and Fig. 5 gives 
the velocity at the photosphere as a function of time and for the 8.38 M© model. In Fig. 4 
only the intrinsically produced ^^Ni was considered, however, the energy release from ^^Ni 
decay has little effect on the expansion velocities, especially for the cases with high explosion 
energy. The fact that speeds as high as log (f /cms~^) = 9.5 were seen on SN 2003dh on day 
10 (Hjorth et al. 2003) illustrates the need for at least some portion of the ejecta to have 
equivalent isotropic energies above 4x10^^ erg. However, most of the models can give the 
milder requirement of log {v/cms~^) ~ 9.0 on day 30. 

4. A "TWO-DIMENSIONAL" MODEL 

Only the most energetic symmetric explosions, E > 10^^ erg, rise rapidly enough, make 
sufficient ^^Ni, and decay rapidly enough to resemble SN 2003dh. The velocities in these 
models are also consistent with what was seen (Fig. 5), but the energies strain credibility and 
even the most energetic model would not make a powerful gamma-ray burst. The correct 
model must be asymmetric. 

Ideally, one would like to explore the coupled GRB-supernova explosion in two-dimensional 
code that couples radiation transport, gamma-ray deposition, and special relativistic hydro- 
dynamics. Such calculations will surely be done, but as an expedient and for clarity in 
exposition, we consider here a simple, composite toy model rendered out of combinations of 
our one-dimensional models. 

We rely here on experience with the collapsar model, especially MacFadyen & Woosley 
(1999) and Zhang, Woosley, & MacFadyen (2002). A black hole forms in the middle of the 
star and, after a few seconds while the polar accretion rate declines, launches relativistic jets 
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along both axes. These jets penetrate the star and ultimately make the GRB, but they do 
not, by themselves, make a supernova. The peak brightness of a Type I supernova of any 
subclass is measured by the ^^Ni that it ejects. Both the jets and the lateral shocks they 
launch make almost none. Lacking additional sources of power, the supernova accompanying 
a collapsar-produced GRB would he weak and nearly invisible. 

That additional source of power is the wind off the accretion disk (MacFadyen & Woosley 
1999; MacFadyen 2002; Pruet, Woosley, & Hoffman 2003). The mass ejected is of the same 
order as the mass accreted by the black hole, i.e., ~ 1 Mq. It's energy is uncertain, but could 
easily be ~ 10^^ erg, that is the binding energy of a solar mass of material at the last stable 
orbit around a Kerr black hole times a few percent. A similar amount of energy is released 
by the reassembly of 0.5 M0 of nucleons from the disk into bound nuclei. Indeed the energy 
in this ''wind" may considerably exceed the energy in the GRB-producing jets themselves. Its 
composition is likely to be chiefly ^^Ni (Pruet, Woosley, & Hoffman 2003) mixed with the 
helium, oxygen, and other heavy elements that make up the star. 

This wind drives a highly asymmetric explosion (MacFadyen & Woosley 1999). In fact, 
to first order, it blows two conical-shaped "wedges" out of the star, each along the rotational 
axis. In the case of SN 2003dh, because we saw the GRB, one of these inverted cones was 
directed straight at us. The equatorial regions are partly ejected, but partly fall into the 
hole, continuing to power the jet for some time after the initial explosion. The opening angle 
of the conical wedges is unknown, but certainly greater than the ~ 5° opening of the GRB 
jet itself, yet probably small enough to contain only a fraction of the stellar mass. Here we 
will use 45° as an example. This implies that 1 — cos 6' = 29 % of the star is ejected at very 
high velocity. Within this 45° we assume a Gaussian distribution of kinetic energies between 
160 X 10^^ erg and 40 x 10^^ erg (Fig. 6) as a function of angle. We further assume a total 
production of 0.55 Mq of ^^Ni, 90% of which comes out in the high velocity (well-mixed) 
wedges, and 10 % of which stays behind in the low velocity ejecta. 

The composite light curves were computed by assuming an explosion energy as a function 
of colatitude, 6: 



This was used to determine the flux as function of 6 by interpolation in the grid of one- 
dimensional light curves. The one-dimensional models were calculated with an effective 




(1) 



where 9i is chosen such that Eiji/A) = 40 xlO 



erg. I.e., 
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mass of ^^Ni equal to 1.5 Mq (last panel of Fig. 3), i.e., a mass such that 29 % of it gave the 
actual ^^Ni mass in the high-velocity ejecta, 0.44 Mq. 

These contributions were then integrated over the projection of the sphere in the direc- 
tion of the observer, along the pole (Fig. 7). Such a procedure does not take into account the 
"shadowing" effect of the fast moving polar ejecta on the slower equatorial ejecta. Owing 
to its fast expansion, the polar ejecta quickly become optically thin and the late-time light 
curve from the equatorial ejecta should not be significantly affected. We also do not take 
into account that the cone of high explosion energy may trap 7-rays less efficiently than a 
full sphere and thus cool down faster, or similarly, that the polar "hole" punched into the 
explosion by the jet could allow photons to escape faster in this direction, or that a larger 
fraction of the low-energy ejecta than corresponds to the projected surface area would be- 
come visible once the polar ejecta are optically thin. These truly two-dimensional effects 
may make the low-E part of the ejecta a little brighter early on and decline a bit faster at 
late times before the spherical model becomes optically thin. 

The composite light curve corresponding to these assumptions is given in Fig. 8. During 
the first month, the luminosity is given entirely by the high- velocity, nickel-rich ejecta. Be- 
cause of the high velocity and complete mixing the rise to maximum is very rapid. Because 
the material becomes thin and gamma-rays from radioactive decay escape, the decline from 
maximum is also fast. The decline of the light curve might be even faster than suggested by 
the composite model, because 7-rays can escape easier from a polar wedge than from a full 
sphere that has been used for modeling the light curves. 

At late times though, the slower moving, nickel-poor ejecta dominate the light curve 
because gamma-ray deposition there remains highly efficient. For situations where the energy 
is due chiefly to radioactivity, this sort of light curve is distinctive signature of asymmetry. 
Models that expand slowly would not be so bright at peak and models that expand rapidly 
would not be so bright at late times. The amount of ^^Ni inferred from the peak brightness 
would be quite different from that inferred on the tail. Such seems to have been the case 
for SN 1998bw where a ^^Ni mass of 0.7 was inferred for the peak (Iwamoto et al. 1998) 
and 0.24 Mq at late times (Nakamura et al. 2001). 

5. DISCUSSION 

The optical light curve of the supernova that accompanies a typical GRB is unique in 
two ways. First, it must be disentangled from the "optical afterglow", which is dominated 
at early times by the relativistic GRB-producing jet interacting with the circum-source 
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medium. There will also be optical emission from slower moving, but still highly energetic 
supernova ejecta running into this medium. Continued, possibly time-variable output from 
the central engine also contributes at a declining rate to the afterglow at all wavelengths 
(Zhang, Woosley, & MacFadyen 2002). Whether one counts these afterglows as "supernova" 
or something else is largely a matter of taste. Supernova shock interactions have often 
been considered part of the light curve (Leibundgut 1994). They are probably the origin of 
variations in the early optical emission seen during the first few days of SN 2003dh (Willingale 
et al. 2003). Even with complete mixing, enormous explosion energy, and a large mass of 
^^Ni, a supernova powered by radioactivity would not reach maximum light in two days. 

Second, the supernova is grossly asymmetric. We have attempted to account for this 
by merging the results of several spherically symmetric models of various energies. Even so, 
large energies are required in the composite model, far above the GRB jet energy estimated 
by Frail et al. (2001). This energy must be provided by another source which, in the coUapsar 
model, is the disk wind. Any attempt to model SN 2003dh without the formation of a black 
hole and disk must find a different way to produce the necessary ^^Ni and energy. The only 
other potential candidate at this time is a millisecond magnetar (Wheeler, Yi, Hoflich, & 
Wang 2000). Provided the putative pulsar can survive accretion in a massive progenitor star 
long enough, it will need to deposit ~ 10^^ erg in ~ 1 s in order to make ~ 0.5 Mq of ^^Ni in 
a high mass progenitor (making this much in a low mass model is impossible). The ^^Ni is 
required by the light curve and the high energy is required to make the star move fast enough 
to explain the rapid post-maximum decline of the supernova. The one-second time scale is 
set by the requirement that the ^^Ni be made by a shock that raises the temperature to 
greater than 5x10^ K as is necessary to produce ^^Ni out of lighter elements. Concentrating 
the pulsar's energy in a smaller solid angle, as has been done here, will give a smaller ^^Ni 
mass no greater than 1.1 Mq times that solid angle (Table 3). In addition, the pulsar would 
need to focus some part of this energy into a very relativistic narrow jet (in order to make 
the GRB) at a time when the accretion rate is very high. 

These potential difficulties lead us to favor the coUapsar model. Based upon this model, 
in which the ^^Ni is provided by the disk wind in an explosion that is highly asymmetric we 
have developed a composite working model in which the total explosion energy is 26 xlO^^ erg 
and the mass of ^^Ni is 0.55 Mq. Further observations to test the predicted light curve (Fig. 8) 
are encouraged. 

This research has been supported by NASA (NAG5-8128, NAG5-12036, and MIT- 
292701) and the DOE Program for Scientific Discovery through Advanced Computing (Sci- 
DAC; DE-FC02-01ER41176). AH has been supported in part by the Department of Energy 
under grant B341495 to the Center for Astrophysical Thermonuclear Flashes at the Univer- 
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Table 1. Properties^ of the 15 Mq progenitor model without mass loss 





m 


r 


J 




(Mo) 


(cm) 


(erg s) 


iron core 


1.95 


2.58x10^ 


1.23x10^° 


Si core 


2.61 


4.99x10^ 


2.35x10^° 


Ne/Mg/0 core 


2.95 


7.06 xlO^ 


2.68x10^° 


C/0 core 


8.56 


5.16x10^ 


2.24x10^1 


star / He core 


15.00 


8.80x10^° 


1.00x10^2 



'^Enclosed mass, m, radius, r, and enclosed angu- 
lar momentum, J, of the progenitor model at core 
collapse or the outer boundaries of the indicated 
cores. 
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Table 2. Properties of the 8.38 Mq progenitor model with mass loss 





m 


r 


J 




(Mo) 


(cm) 


(erg s) 


iron core 


1.46 


1.63x10^ 


4.20x10^9 


Si core 


1.77 


5.06x10^ 


5.57x10^9 


Ne/Mg/0 core 


1.88 


6.77x10^ 


5.96x10^9 


C/0 core 


4.71 


6.12x109 


3.50x10^° 


star / He core 


8.38 


8.09x10^° 


1.38x10^1 
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Table 3. Parameters of the Explosions 



Mass (Mq) 


8.38 


15 


15 


Piston (M0) 


1.45 


1.93 


2.75 


KE (10^1 erg) 




5^Ni (Mq) 




1.25 


0.158 


0.040^^ 


0.024^ 


2.5 


0.193 


0.614^^ 


0.225 


5 


0.229 


0.707 


0.254 


10 


0.244 


0.815 


0.303 


20 


0.230 


0.938 


0.377 


40 


0.189 


1.045 


0.487 


80 


0.168 


1.083 


0.583 


160 


0.170 


0.955 


0.625 



'^The ^^Ni mass is reduced by fallback 
after mild mixing. For the piston loca- 
tion at 1.93 M0 the ^^Ni masses before fall- 
back are 0.558 Mq and 0.623 Mq; the to- 
tal fallback mass is 6.889 Mq and 0.050 M© 
(mass coordinates 8.819 Mq and 1.980 Mq). 
For the piston location at 2.75 Mq the ^^Ni 
mass without fallback is 0.192 Mq, the fall- 
back mass is 5.580 Mq (mass coordinate 
8.330 Mq). 
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Fig. 1. — Left: Composition of the 8.38 M© star with 1.25 xlO erg explosion energy 100s 
after core collapse, but prior to any mixing. Right: Same plot, but for the 160 x 10^^ erg 
explosion. Note that 0.8 Mq has been removed from the surface of this calculation because 
it exceeded 10^°cms~^ and posed problems for the non-relativistic hydro-code. Also note 
the large mass of photodisintegrated matter (^He) in the inner regions near the mass cut. 
This limits the production of ^^Ni in the model. 



- 16 - 




5 10 15 20 25 

time / days 

Fig. 2. — Light curves for "moderately" mixed models derived from the high mass progenitor 
(15 Mq) and the large piston mass (2.75 M0). The explosions were calculated in spherical 
symmetry with kinetic energies at infinity of 160, 80, 40, 20, 10, 5, 2.5, and 1.25 xlO^^ erg and 
the light curves include just the mass of ^^Ni produced explosively in the model (Table 3). 
Note the pre-maximum non-monotonic evolution of the most energetic models resulting from 
helium and heavy element recombination. For the most energetic explosions the initial rise 
for the first few hours is a numerical artifact of removing zones moving faster than 1/3 c. 
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Fig. 3. — Light curves from the models for the low mass progenitor (8.4 Mq). Calculated 
in spherical symmetry with kinetic energies at infinity of 160, 80, 40, 20, 10, 5, 2.5, and 
1.25 xlO^^ erg. See also Fig. 2. Top Left: Moderately mixed with the ^^Ni yields in Table 3. 
Top Right: Same models as top left, but completely mixed. Bottom Left: Completely mixed 
models with a normalized abundance of ^^Ni of 0.5 M© in all cases. The peaks of the light 
curves are brighter owing to the increased mass of ^^Ni and the pre-maximum variability is 
less apparent. Bottom Right: Same mixed calculations normalized to a ^^Ni mass of 1.5 Mq, 
as is appropriate for the high-velocity ejecta. 
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Fig. 4. — Left: Final velocity as a function of mass for a piston location of 1.93 M0 in the 
15 M0 model. Moderate mixing and only the ^^Ni produced intrinsically by the explosion 
were used. Note fallback of ~ 9 Mq for the lowest explosion energy. Right: Same plot, but 
for the 8.38 Mo model 



- 19 - 




u 8.5 - 
_o 

cu 
> 



8.0 f . . . . I . . 1 . I . . I. . I . .1 . . I . . . I. I . . . . 1 f . . . . I . . . . I . .1 . . I . . .1 . I . . . .1 1 . . . . ' 
5 10 15 20 25 5 1 1 5 20 25 

time / days time / days 

Fig. 5.— Left: Velocity at the photosphere for the 160, 80, 40, 20, 10, 5, 2.5, and 
1.25 X 10^^ erg explosions of the 8.38 star with moderate mixing and only its intrinsic 
^^Ni production. Note that the velocity drops to zero and becomes undefined thereafter 
when the photosphere reaches the center of the ejecta. (A line of sight through the center 
would have an optical thickness of 4/3 at this point.) Also note that zones moving faster 
than lO^'' cms~^ have been cut off from the calculations; this affects the velocity at the pho- 
tosphere for about half a day for the most energetic explosion, and increasingly earlier times 
for lower energies. Right: Same plot, but for ejecta whose ^^Ni mass has been increased 
to 1.5 Mq and the ejecta were "thoroughly mixed". Extra heating by the increased ^^Ni 
abundance at large radii maintains an extended photosphere longer. 
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Fig. 6. — Contributions to the composite model. The figure gives the equivalent isotropic 
energy along a given angle as a function of angle. For angles less than 45° the energy varies 
from 160 xlO^^ erg to 40 xlO^^ erg and models normalized to 1.5 M0 of ^^Ni were used, being 
"thorough" mixed. For angles greater than 45° the explosion energy is 1.25 x 10^^ erg weak; 
only the intrinsic ^^Ni production (Table 2) and mild mixing was employed. 
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Fig. 7. — Decomposition of the light curve. Assuming, a (bipolar) jet of opening angle 29 
(e.g., 6 = 45° as depicted above), each polar cap comprises a solid angle of 27r (1 — cos6'), 
i.e., the bipolar jet flows out at a fraction 1 — cos 6* of the total solid angle of Air. This 
means 1 — cos^ of the isotropic explosion contributes to such an explosion and its energy. 
However, an observer looking along the axis of the jet sees the projected surface area vr (sin 9)'^ 
out of 271 for the full circle, the full projected sphere, i.e., a fraction sin^ ^ of the isotropic 
light curves contributes to apparent light curve for a such positioned observer. Note that 
for 0° < ^ < 90°: sin^ ^ > 1 — cos 5, i.e., the observer along the axis is favored to see an 
apparently brighter light curve than would correspond to the total energy of a composed light 
curve (assuming the explosion is stronger in polar direction than in equatorial direction). 
This is because material at high colatitude, 9, appears under a high inclination angle. 
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Fig. 8. — Composite light curve resulting from a combination of well-mixed models (Fig. 3) 
having a distribution of kinetic energies as given in Fig. 6. For polar angles less than 45° 
the three highest energy models in Fig. 3 were employed which each have an effective ^^Ni 
mass of 1.5 M0. Together, with appropriate angular weighting (see text), these give the 
curve labeled "high-E". For angles greater than 45° the model with energy 1.25 x 10^^ erg in 
Fig. 3 was used. The ^^Ni mass in this component was left at 0.15 Mq. This curve, with 
appropriate weighting, is labeled "low-E" . The dark solid line is the total. The total energy 
in the composite is 26 xlO^^ erg and the total mass of ^^Ni is 0.55 Mq. 



